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Abstract 

Q_i' We reanalyze the supersymmetric leptogenesis model by one of the 

authors (HM) and Yanagida based on L = H u flat direction in detail. 
D , We point out that the appropriate amount of baryon asymmetry can 

be generated in this model with the neutrino mass matrix consistent 
with the atmospheric neutrino oscillation and solutions to the solar 
neutrino problem, preferably the small angle MSW solution. The 
reheating temperature can be low enough to avoid the cosmological 
gravitino problem. This model is the minimal one because it does 
not rely on any new physics beyond supersymmetry and Majorana 
neutrino masses. 



* This work was supported in part by the U.S. Department of Energy under Con- 
tracts DE-AC03-76SF00098, in part by the National Science Foundation under grants 
PHY-95-13835 and PHY-95-14797. TM was also supported by the Marvin L. Goldberger 
Membership. HM was also supported by Alfred P. Sloan Foundation. 



1. Introduction 

Cosmic baryon asymmetry has been one of the biggest puzzles in mod- 
ern cosmology. It is hoped to understand it in terms of microscopic physics 
and evolution of Universe, if the well-known Sakharov's three conditions are 
satisfied: (1) existence of baryon-number violation, (2) existence of CP vio- 
lation, and (3) departure from thermal equilibrium. A stringent constraint 
on the baryogenesis models is posed by the apparent lack of baryon-number 
violation. Recent data from Soudan-II and SuperKamiokande experiments 
|lj made the constraints on proton and neutron decay yet more stronger. 
Especially given the constraint on p — > e + 7r°, which requires the mass of the 
X boson in S77(5) grand unified theories to be above 10 15 GeV or so, it is be- 
coming increasingly difficult to construct models of baryogenesis with explicit 
baryon-number violation. On the other hand, the recent evidence for oscilla- 
tion in atmospheric neutrinos from SuperKamiokande M, Soudan-II [fj] and 
MACRO experiments strongly suggest small but finite neutrino masses 
with Am 2 = 10~ 3 -10 -2 eV 2 . Such small mass scales are unlikely to arise for 
Dirac neutrinos,[] while are naturally understood in terms of so-called seesaw 
mechanism [||] which in turn requires Majorana neutrinos. Therefore, given 
the theoretical prejudice, there is a strong evidence for lepton-number viola- 
tion at the mass scale of the right-handed neutrinos Mr ~ 10 14 -10 15 GeV. 

Fukugita and Yanagida J7| suggested that the lepton-number violation 
may be enough for baryogenesis because the electroweak anomaly violates 
B + L H and, once lepton number is violated, so is baryon number. If a 
finite lepton asymmetry is created in the Early Universe, it can be partially 
converted to baryon asymmetry without a need for an explicit baryon-number 
violating process which could lead to too-rapid proton decay. This scenario, 
called leptogenesis, is hence probably the most attractive possibility of baryo- 
genesis in the current situations. The thermally produced right-handed neu- 
trinos decay out of equilibrium and the CP-violation in the neutrino Yukawa 
matrix causes a difference between the decay rates into leptons and anti- 
leptons in this model. The lightest right-handed neutrino is typically as 
heavy as 10 10 GeV in many models and therefore the reheating temperature 
after inflation needs to be higher than that (see || for a recent discussion). 

On the other hand, the large hierarchy between the electroweak scale 
and the scale of lepton number violation Mr would require a stabilization 

TAn alternative possibility was pointed out recently using large extra dimensions with 
right-handed neutrinos in the bulk [pi. 
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mechanism for the hierarchy. Supersymmetry is the best solution to this 
problem currently available. Therefore it is natural and timely to ask what 
the minimal scheme is for the supersymmetric leptogenesis. 

Supersymmetry, however, causes a cosmological problem when the re- 
heating temperature after the primordial inflation is too high. The gravitinos 
can be produced thermally whose abundance is approximately proportional 
to the reheating temperature. For a typical mass range m 3 / 2 ~ 500 GeV 
favored to stabilize the hierarchy^ the produced gravitinos decay after the 
Big-Bang Nucleosynthesis and they dissociate light elements. To keep the 
theory and observation consistent, there is an upper bound on the gravitino 
abundance, and hence on the reheating temperature of about T RH < 10 9 GeV 
for this mass range (for the latest analysis, see This constraint makes 

it somewhat difficult to produce right-handed neutrino states as required in 
the original model |7|].fj] 

In this letter, we study the minimal supersymmetric leptogenesis. Our 
definition of minimal is the particle content of the Minimal Supersymmet- 
ric Standard Model together with suggested neutrino masses as effective 
dimension-five operators. This way, we avoid reliance on explicit grand- 
unified models or seesaw models. The only dependence on physics beyond 
the standard model is supersymmetry and neutrino masses. We show that 
leptogenesis is possible a la Affleck-Dine [[H] within this minimal framework 



as originally suggested by one of the authors (HM) and Yanagida [15| , within 
the current experimental situations. The reheating temperature can be low 
enough to avoid the gravitino problem because the model does not require 
the production of right-handed neutrinos. 
2. The Model 

We start with the effective dimension-five operators in the superpotential 

Tfl 1 

W = ^|(LA)(LA) = — {UH^L^) (1) 

which are necessary to incorporate finite neutrino masses (and hence lepton- 
number violation) to the Minimal Supersymmetric Standard Model. The 
vacuum expectation value for H u is v u = (H u ) = (174 GeV) sin/5 and we 

* Anomaly mediation of supersymmetry breaking [ |l0| may allow much higher gravitino 
mass keeping most superparticles light, and hence a higher reheating temperature. 

§One can get around this problem, however, if the right-handed neutrinos are produced 
non-thermally, i.e., by the decay of the inflaton [O] or preheating p3|. 
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assume a moderately large tan (3 <; 5 throughout the paper .[] The index i runs 
over three mass eigenstates of neutrinos. It will turn out that the relevant 
direction is most likely the lightest neutrino as we will see later and henceforth 
suppress the index i. Along the flat direction L = H u = (fi/V^, (the factor 
\/2 is necessary to ensure the canonical kinetic term: L*L + H*H U = </>*</>) 
the potential is given by 

Here, the first two terms are from SUSY breaking and we assume m ~ 
A ~ 100 GeV — 1 TeV, while the last term is from superpotential given in 
Eq. ([!]). The mass scale M does not necessarily be the mass scale of right- 
handed neutrinos; it can actually be much higher. Our philosophy is not 
to discuss the precise origin of such operator, but rather use this effective 
operator only. 
3. Initial Amplitude 

Since this direction is D-flat and also F-flat in the vn v — > limit, it can 
have a large amplitude at the end of the primordial inflation. It turns out 
that the resulting baryon asymmetry does not depend on the precise value of 
the initial amplitude. All we need is that the initial amplitude is larger than 
\fmM as we will see later. We would like to discuss various mechanisms to 
generate a large amplitude proposed in the literature. We have the expansion 
rate of the Universe during the primordial inflation of if ,„/ ~ 10 12 -10 14 GeV 
as suggested by popular models |16| . 



To discuss the amplitude in the Early Universe, one needs to pay special 
attention to the fact that the expanding Universe itself breaks global super- 
symmetry and may modify the scalar potential (TTS|. The explicit form of 



the supersymmetry breaking, however, depends on the details of the Kahler 
potential and hence is model dependent. Therefore we consider three possi- 
bilities in this letter. (1) No modification from the flat-space potential, which 
occurs at the tree- level of the no-scale type supergravity fl9fl , (2) loop-level 



modifications from the flat-space potential JI5 |, and (3) tree-level modifica- 
tions from the flat-space potential [0 . In all cases, we find it is easy to have 
an initial amplitude necessary for the model: 4>i n j ^> \JmM as we will see 
later. 



^For smaller tan (3 ^ 1.4 allowed by the perturbative top Yukawa coupling, the resulting 
lepton asymmetry changes only up to a factor of two. 
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In case (1), one idea to generate a large initial amplitude is to start with a 
"chaotic" initial condition (T^],| where all scalar fields start with large values. 
During the course of inflation, most of the fields are diluted exponentially due 
to their potential, but some remain if their potential is flat enough. Along 
our flat direction, the potential is eventually dominated by |0| 6 term once 
the amplitude becomes small enough whatever the physics at high energy 
scale 3> M is. The evolution is given by the slow-rolling approximation, and 
we can solve the equation of motion 

3#i B /0+ -j^MV = 0. (3) 

For a more-or-less constant expansion rate H = Hi n j during the inflation, its 
solution is 

1 I 1 

WW) ~ WW) ~ M*Hw ' 
Therefore with an e-folding of iV = H in ft, 

mt) = \<P\-'(o) + N/(Mmi f y (5) 

Clearly, the amplitude at the end of inflation is almost predicted as long as 
|0(O)| is large enough, 

M 2 H 2 

\mt) = (6) 

The e-folding N required to solve the horizon and flatness problems is iV > 
100 [0. 

The cases (2) and (3) use a possible negative mass squared of order 
-H 2 nf /16n 2 U or -H 2 nf @. Minimizing 

V = -cHl f \^ + (7) 

with c ~ 0(1) or 0(1/16tt 2 ), we find -cH 2 nf( j)+ (3/4M 2 )|0| 4 = 0, or |0| 4 = 
4 C ^M 2 /3. 

Both above arguments are classical. One may worry that the quan- 
tum fluctuation during the de Sitter expansion of Universe may modify the 
estimate. An effective Hawking-Gibbons temperature is given by Thg = 



"We do not need V ~ A/p lanck a s suggested in [n7|, however. 
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H in f/2iT and the fields acquire expectation values which correspond to the 
potential energy density of the order of the Tfj G pOf . This effect suggests 



M 2 Hf n j, which is typically smaller than the other estimates for the 
range of M of our interest (see below). This smallness of the quantum fluc- 
tuation therefore justifies the classical treatment above for the flat direction 
because the fluctuation around the classical value can be neglected. 

Based on the above discussions, we take an order of magnitude estimate 
of the amplitude of at the end of the inflation ti n f. 

<t>t nf = \<f>(t mf )\ A ~ M 2 H 2 nf , (8) 

or slightly smaller for the rest of the analysis. The result turns out not to 
depend on the precise value of the initial amplitude; the only requirement is 
that (pi n f ^> v mM. 
4. Evolution 

Given the above initial condition, we follow the cosmological evolution 
after the end of the inflation. While the inflaton oscillates around the min- 
imum of its potential, the flat direction gets diluted. For case (1), the 
equation of motion of the when the potential is dominated by 6 term is 
given by 

<») 

assuming the matter dominated expansion H = 2/31 For case (2), there 
is an additional term of order — 16 ^ 2f 2 in the l.h.s. of the equation, and 
similarly — 4-0 for case (3). Despite the apparent difference between three 
cases, turns out to track yMH until the expansion rate slows down to 
H ~ m as we will explain below. 

We first discuss case (1) with flat-space potential because of its simplicity. 
The behavior of |0| depends on which term dominates its potential, since the 
potential changes its shape at 



|0| ^V^M = (p LG - (10) 

At the time when |0| ~ <Plg-, the lepton number gets fixed (we refer to this 
point as LG, for LeptoGenesis). For |0| > iG , |0| 6 term dominates the 
potential and the equation of motion (^) should be used. Using the slow-roll 
approximation |0| -C |20/i|, the solution is given by 



t 2 - 1 2 

,-4 , L inf 

inf AM 2 
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and it is easy to check that the motion is critically damped when (f) in f < 
y 1 MH. Therefore, independent of the precise initial value of 4>i n f, the mo- 
tion starts only when <pi n f ~ y MH. After (f) starts to move, one can use two 
extreme approximations to study the dilution of <fi: the slow-roll approxima- 
tion and the virial motion. With the slow-roll approximation, the solution 
(|TT| ) above dilutes as 4 cx R~ 3 oc H 2 . The opposite limit of the virial (os- 
cillatory) motion with 6 potential gives <j) 2 /2 = 3V. Then the equation of 
motion + 3H<p + V = can be rewritten in terms of the energy density 
P(j> = 2 /2 + V as = —3H(p 2 = —^Hpcf,, and hence p^ oc R~ 9 / 2 . Therefore 
we again obtain 4 oc R~ 3 oc H 2 . The true motion is somewhere between the 
slow-roll and virial limits but the agreement of both extremes suggests that 
this is the true power law of dilution. Then the amplitude always traces 
y 1 MH. On the other hand, once \<p\ ^ 4>lg, (ft has a quadratic potential. In 
this case, the particle picture is valid for and (f> dilutes as \<fi\ 2 oc R~ 3 . 

For a low reheating temperature (Trh ~ 10 10 GeV), \<p\ ~ (pic is realized 
in the inflaton-dominated Universe. Because 4> tracks y 1 MH, this occurs 
when H ~ Hlg — m - This corresponds to the inflaton energy density 

p^t LG )~3m 2 Ml (12) 

where M* = M Planck / 'yf&K = 2.4 x 10 18 GeV is the reduced Planck mass. 

Lepton number density at this time is easily estimated. Define the (f>- 
number density = i(<f)*<j) — <f)*<f)), which is related to the lepton number 
density as til = \ n <f>, then its equation of motion is given by 

n + 3Hn^ = i ( - h.c] = ^M^ 4 )- (13) 



V d(j) J M 



As a result, we find that 



|(it%) = i^lm(0 4 ), (14) 



and hence 

r°° A 

{R Z n^){t LG ) = j o dtR 3 —lm(^). (15) 

For t < t^c in the integral, 6 dominates the potential and 4 oc R~ 3 . The 
integral for this time range just gives the time interval tzG. For t > t^Q, 
(p 2 dominates the potential and 4 oc R~ 6 . The integral for this time range 
gives the same contribution as the previous one. Therefore the resulting 
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lepton asymmetry is dominated by the late time contribution at H ~ m and 
rather insensitive to the history before then. For matter-dominated Universe, 
tj_,Q — ^H^q, and hence 

44 

MtLc) - (16) 

By substituting 4 ~ m 2 M 2 e At6 and A ~ m, 

4 

n^(t LG ) ~ -m 2 M sin40. (17) 

The size of the resulting lepton number is determined by the initial phase of 
the field 9. This is because the potential has "valleys" along 9 = ir/4, 37r/4, 57r/4, 77r/4 
and "ridges" along = 0, tc/2, it, 3tc/2 and the initial phase generates a rota- 
tional motion of on the complex plane when it rolls down the "slope" from 
the "ridge." In this sense, the lepton number is determined by a "sponta- 
neous" CP violation because of the random phase the field acquired during 
the inflation. 

When the inflaton decays, the Universe is reheated. With the instanta- 
neous decay approximation, p^itRu) = ^q9*Tr H , and we obtain 



Rrh\ 3 _ Pf(t LG ) 3m 2 M, 2 



(18) 



RlgJ P^rh) ^9*Tr H 
Therefore, the number density at the reheating is given by 

rictus) = n^t LG ) (^) ~ g^ sin4 ^ x 1)^ T ^' (19) 

or using s{t RH ) = ^-g*T^ H , we obtain 

MT RH . 

The above estimate is valid as long as |0| ~ 4>lg is realized in the inflaton- 
dominated Universe. This is the case if 3m 2 M 2 > ^q9*Trh, or equivalently, 
Trh ~ 10 10 GeV for m ~ 100 GeV. For higher reheating temperature, how- 
ever, |0| ~ <Plg happens after the reheating. At the end of the inflation, H = 
H in f, and the energy density of the oscillating inflaton ip is p^ = 3Hf n jM 2 . 
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Using p^t RH ) = ^g*T RH , we find {R inf /R RH f ~ (^T RH ) /(3H 2 nf M 2 ). At 
the same time, the flat direction is diluted as 



nt RH )~M*Hi f (^-\ ^m*^ (21) 



' mf \R RH ) ~ 3M2 " 30' 

which is still larger than for > 10 10 GeV. 

At the time when |0| 4 ~ m 2 M 2 , the lepton number gets fixed. The scale 
factor at this point is determined by 

= (M ~ (22) 

\Rlg/ \<PrhJ ^9*Tr H 
The radiation energy density at t — t^G is therefore given by 

A-(t M ) = { To S.T RH ) {-^) ~ 3m M, (^-j . (23) 



This gives the expansion rate at this point 

3m 2 M 2 



1/6 



HLG~m[^—t) ■ (24) 

The evolution of n^, is solved as in the previous case. For radiation-dominated 
Universe, tic = \Hlg-, and hence Eq. (|15D gives 




n^tLo) ~ -m^M a * sin 40, (25) 



30 



where we used 4 ~ m 2 M 2 e ild and v4 ~ m. If we naively adopt the above 
relation, is bigger than m 2 M ~ ^I^lgI 2 f° r @ ~ ^(1)- However, at this 
point, the motion allows a particle interpretation because the potential 
is dominated by the quadratic piece, and the maximum possible number of 
particles is given by ~ m|0| 2 . If we take the above estimate literally, 
may be bigger than m|0| 2 : a contradiction. Then the conclusion is that 
has already saturated before this stage. Therefore, some of the previous 
discussions for T R h ~ (mM,) 1 ' 2 are somewhat irrelevant, and we can simply 
estimate at ~ <Plg to be ~ m 2 M. Then, using the entropy at tiG- 

, x / x fRRH\ 3 4m 2 M 2 , . 

s(t LG ) = s(t RH ) ~^=-^, (26) 

V t^LG / 1 RH 
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we obtain a similar 0-to-entropy ratio as given in Eq. (|2H| ) up to an 0(1) 
coefficient. We use the result given in Eq. ( p0|) for the numerical estimation 
below. 

With the estimated n^/s (0), the lepton-to-entropy ratio is given by 

v - Ul n $ lMT RH ■ AQ lvlT RII . , , 

Y L = — = — ~ — — — sin 46*= — — sin 40. 27 

Since the lepton number is inversely proportional to the mass of the neutrino, 
it is dominated by the lightest mass eigenvalue m Ul , possibly mostly the 
electron neutrino state. 

So far the discussion used only the flat-space potential (case (1)). For 
case (2), the discussion is modified only slightly. The field amplitude dilutes 
as (ft ~ v MH because the negative mass squared —H 2 /16it 2 is smaller than 
the expansion rate and hence does not modify the discussion with flat-space 
potential. The main issue here is if the phase direction of the field gets 
settled to the "valley." Because of loop-suppressed effect of the expansion 
in the potential, the potential along the phase direction is given by ((A + 
cif/167r 2 )(/> 4 + c.c.)/8M with c = 0(1), and the curvature along the phase 
direction is of the order of cH<p 2 /167r 2 M ~ cH 2 /16tt 2 . This is always smaller 
than the expansion rate and hence the motion along the phase direction is 
critically damped. Therefore the initial phase of the potential is kept intact 
until the expansion slows down to H ~ m when the flat-space potential is 
recovered and the field starts rolling down the slope from the "ridge." The 
estimate of the resulting lepton asymmetry is hence the same as in the case 
(1) with the flat-space potential. 

The discussion is different for case (3) where the effect of the expansion 
introduces a large supersymmetry breaking effect of 0(H) in the potential. 
This case was discussed in detail in ||21|| , and we summarize the situation 
briefly here. The field amplitude always tracks the minimum of the potential 
~ \/MH because the curvature — H 2 \(p\ 2 is of the same order of magnitude 
as the expansion rate itself and field always rolls down to the minimum at any 
given time. Therefore the field amplitude is the same as other two cases even 
though the reason for it is quite different. Similarly, the potential along the 
phase direction ((A+cif )0 4 + c.c.)/8M with c = 0(1) produces the curvature 
along the phase direction is also of the order of cH<p 2 /M ~ H 2 and hence the 
field settles to the "valley" and tracks it. However, the potential can generate 
a rotational motion when H becomes smaller than A because of a possible 
relative phase between two terms A and cH. The phase for the "valleys" 
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quickly shifts to that of the flat-space potential at this moment and the field 
which tracked the valley with the H effect starts moving to the new valley 
with the A term only. The estimate of the rotational motion is basically the 
same as the earlier two cases and hence the resulting lepton asymmetry as 
well. What is different from the previous cases is that the CP violation is 
in the relative phase between A and cH and originates from the microscopic 
physics rather than a spontaneous random phase. 

To summarize, in all cases (1), (2), (3) considered, the final lepton asym- 
metry from the motion of the .D-flat direction L = H u is given by Eq. fl27|) 
which confirms the original estimate in fI5 |. Therefore the result is quite 
robust and model independent. 
5. Estimate of Baryon Asymmetry 

In order to proceed to a numerical estimate of the generated lepton asym- 
metry, we need to have some idea on the lightest neutrino mass from the 
currently available data. Here we have to rely on certain assumptions. 

We assume hierarchical neutrino masses similar to quark and charged 
lepton masses, and hence that Am 2 from neutrino oscillation give larger of 
the mass eigenvalues between two neutrino states relevant to the particular 
oscillation mode. The atmospheric neutrino data require that the largest 
two eigenvalues have Am 2 = 10 _3 -10~ 2 eV 2 , and hence according to our 
assumption, m U3 = 0.03-0.1 eV. The small angle MSW solution suggests 
Am 2 ~ 4-10 x 10~ 6 eV 2 , which should give the second lightest mass eigen- 
value m U2 ~ 0.0020-0.0032 eV and the mixing angle sin 2 29 MNS ~lx 10~ 3 - 



1 x 10 [22]. The estimate of the lightest mass eigenvalue depends on the 



assumptions on the mass matrix. If the MNS (Maki-Nakagawa-Sakata) mix- 
ing angle |23| is solely due to the neutrino Majorana matrix, the lightest 
eigenvalue is likely to be around m ui ~ m vfi\iNS ~ (0.5-8) x 10~ 6 eV. 
This approximate relation between the mass eigenvalue and the angle is 
quite generic unless most of the mixing is attributed to the charged lep- 
ton sector. For instance the models based on coset-space family unification 
on E 7 /SU(5) x U(l) 3 [f24}[ , string-inspired anomalous U(l) f2"5|] , and SU(5) 



grand unified model with Abelian flavor symmetry [^B| all give this approxi- 
mate relation. However if the mixing angle is coming from the charged lepton 
sector, the mass eigenvalue can be smaller. 

For the large angle MSW solution, m U2 ~ 0.003-0.01 eV, and the lightest 
eigenvalue is likely to be not much smaller than m V2 if the large MNS angle 
comes from the neutrino mass matrix, but can be much smaller if the MNS 
angle originates in the charged lepton mass matrix. This might be considered 
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a fine-tuning however because the determinant of the mass matrix (either 
charged lepton or neutrino) should be much smaller than the typical size of 
the elements. The vacuum solution most likely gives m Vl ~ m U2 ~ 10~ 5 eV, 
which may also need a fine-tuning in the mass matrices. For the purpose of 
the estimation below, we take m Ul ~ 2.5 x 10~ 6 eV, preferred by the small 
angle MSW solution, but the result can be easily scaled according to the 
value of m vx . 

The reheating temperature cannot be arbitrarily high because of the cos- 
mological constraints on the gravitino production. For the range of gravitino 
mass usually considered, m 3 / 2 — 500 GeV or so, the primordial gravitinos 
will decay with the lifetime of r 3 / 2 — 2 x 10 5 sec and therefore can poten- 
tially destroy the light elements synthesized by Big-Bang Nucleosynthesis. 
The most robust constraint is from the thermally produced gravitino, which 
requires the reheating temperature to be less than Trh ^ 10 9 GeV ITT 



Another constraint is obtained from the non-thermal production after the 
inflationary epoch, which may give a more stringent constraint on the re- 



heating temperature |[27|| . However, the number density of the non-thermally 
produced gravitino is model-dependent. For example, for the chaotic inflation 
model with the inflaton superpotential of W in fl aton ~ lM inflaton <j) 2 inflaton , the 
gravitino to entropy ratio is of order M in fl aton TRH I 'M% , which is smaller than 
that of the thermally produced gravitino. Because of the model-dependence, 
we do not include the constraint from the non-thermally produced gravitino 
into our discussion, and we take the canonical value of Trh = 3 x 10 8 GeV 
for the estimation below. 



Now come the numerical estimates. We find the lepton asymmetry (|27|) 
to be 

Y L ~ 1.1 x 10-° ( 25 * 10 "° CV ) ( r " ,W (28) 



3 x 10 8 GeV 



From the chemical equilibrium between lepton and baryon numbers, we ob- 



tain Yb = 0.351/, 28 , while the number density of photons now is related 



**Moduli fields may also be dangerous. In particular, if they have large initial amplitudes 
of (3(M*), their primordial number density becomes so large that the Big-Bang Nucleosyn- 
thesis is seriously damaged. The initial amplitude may be suppressed, for example, if the 
minimum of the moduli potential is given by a symmetry enhanced point [fl8| . Non-thermal 
production of the moduli may be also important, but is also model-dependent j27j. In par- 
ticular, if the modulus have a large effective mass of 0(Hi n f) during the inflation, or if 
the modulus is conformally coupled, then its non-thermal production after the inflation is 
not effective. 
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to the entropy density by n 7 /s = (410.89 cm 3 )/(2892.4 cm 3 ), and in the 
end we find 



7/ 



— ~ 2.46Y L ~ 2.6 x 10 



n. 



-10 



( 



2.5 x 10" 6 eV 



m 



) 



( 



3 x 10 8 GeV 



Trh 



) 



sin AO. 



(29) 



Recall that we used the preferred range by the small angle MSW solution 
m Ul ~ (0.5-8) x 10~ 6 eV. This estimate should be compared to the value de- 
termined from the Big-Bang Nucleosynthesis r] ~ (4 — 6) x 10 -10 . Therefore 
this mechanism can generate the cosmic baryon asymmetry of the correct 
order of magnitude without conflicting the constraints on the cosmological 
gravitino production within the range of neutrino mass preferred by the atmo- 
spheric neutrino data and the small angle MSW solution to the solar neutrino 
problem. However, one should note that such a small neutrino mass can be 
accommodated even with other solutions to the solar neutrino problem (i.e., 
large angle MSW and vacuum oscillation solutions) by attributing the MNS 
angle to the charged lepton sector or by a fine tuning in the neutrino mass 
matrix. 

6. Conclusion 

To conclude, we have shown that the minimal supersymmetric standard 
model together with Majorana mass of neutrinos can generate the cosmic 
baryon asymmetry without any additional new physics, hence the minimal 
supersymmetric leptogenesis. The required neutrino mass is consistent with 
the mass matrices proposed to explain the atmospheric as well as solar neu- 
trino data with the small angle MSW solution. Other solutions to the solar 
neutrino problem can also be accommodated if the lightest neutrino mass 
eigenvalue is of the order of 10~ 6 eV. The reheating temperature can be low 
enough to avoid cosmological gravitino problem. 
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